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ABSTRACT 


This paper describes a study of low latitude 
middle atmosphere ionization. With data obtained from 
three blunt conductivity probes and one Gerdien con- 
denser, all rocket launched and deployed on parachutes 
during the 1975 ”Antarqui" in PerU; an investigation has 
been conducted into the effects of various ionization 
sources in the 40 to 65 Km altitude range. A comparison 
of the two instrument types suggests that they yield 
consistent results. An observed enhancement of positive 
ion conductivity taking place during the night can be 
explained by an atmospheric effect, with cosmic rays 
being the only source of ionization, only if the ion- 
ion recombination coefficient (a^) is small (<10”^cm^s"^) 
and varies greatly with altitude. More generally accepted 
values of ct,* x 10 cm s ) require an additional 
source of ionization peaking at about 65 Km, and cor- 
responding approximately to the integrated effect of an 
X-ray flux measured on a rocket flown in conjunct on with 
the ionization measurements. I’he physically reasonable 
assumption of an which does not vary with altitude in 
the 50-70 JCm range implies an even greater value and 
a more intense and harder X-ray spectrum. 


vi 



CHAPTER I 
INTRODUCTION 

1 . 1 Middle Atmosphere Ionisation 

Our knowledge of the middle atmosphere (30-80 Km 
altitude range), including the lower ionosphere, is incom- 
plete. This region is characterized by large neutral species 
density and a weakly ionized plasma, yet the fundamental but 
complex mechanisms and atmospheric constituents in this 
region are still only partially understood. Species con- 
centrations and chemical reaction rates, for example, are 
usually given as approximate values and, in many cases, 
only order of magnitude estimates. The uncertainty in these 
atmospheric parameters make it extremely difficult to dis- 
tinguish the relative importance of various atmospheric 
processes . 

Despite the aforementioned difficulties, some 
general observations about ionization sources in the lower 
ionosphere can be made. The ionization of the middle 
atmosphere is produced by several energetic, therefore deeply 
penetrating, energy sources. According to accepted views 
(IVhitten and Poppoff, 1969), the most significant ionization 
sources under undisturbed or "quiet" solar and geophysical 
conditions is the (Solar) Lyman a (1216 A°) ionization of 
NO (Nicolet and Aikin, 1960). If disturbed solar conditions 
exist, the enhanced X-ray flux from the sun becomes important. 
Finally, at lower altitudes during the day and at all altitudes 
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at night galactic cosmic rays which will ionize; all atmo-' 
spheric constituents become predominant (Velinov, 1968). 

The crossover between cosmic ray and Lyman a production 
is usually at 63.5 Km (Hale, 1973). 

Recently, however, a number of investigators 
suggest that an additional nighttime source of ionization 
must be considered: galactic X-rays (Ananthahrishnan, 1969; 

Chilton, 1972). Most of these observations concerned the 
change of VLF signal phase between two monitoring stations. 
It was reported that increased electron production rates 
due to the transit of strong X-ray sources Sco X-1, Tau X-1, 
and the galactic center produced these observed changes in 
VLF phase (Sharma, 1972; Mitra, 1972). This matter is 
still an issue of controversy since other investigators 
contend the increased electron production due to celestial 
X-rays is only on the order of a few percent (Poppoff and 
Whitten, 1969). 

Goldberg (1975) has proposed that some X-ray 
effects may not be due to celestial sources, but Bremsstrah- 
lung produced in the upper ionosphere (F-Region) . 

1 , 2 Importance and Scope of This Research 

This research will contribute to the knowledge of 
the middle atmosphere in a number of ways. Firstly, very 
few studies of this height range have been done at equa- 
torial latitudes . The data employed in this work v;ere 
obtained during the ^Antarqui" D-Region campaign in Peru, 
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in May of 1975. CThis was intended to be an intensive 
investigation of the constituents and processes of this 
region.) This research should then contribute appreciably 
to the understanding of middle atmosphere processes in 
these latitudes. Further, conductivity and density 
data were obtained by two different types of direct sensing 
devices constructed here at The Ionosphere Research Labora- 
tory under the guidance of Dr. L. C. Hale: the Blunt Probe 

and the Gerdien Condenser, both are rocket launched and 
parachute borne instruments. This will only be the second 
time that the conductivity measurements obtained by 
these two experiments, at low latitudes, will be compared. 
(Meteorological balloons or electromagnetic wave techniques 
are not useful in this region.) 

Secondly, these data, combined with reasonable 
assumptions, will be used in solving a simple production- 
loss (mechanism) continuity equation. Since cosmic ray 
production is well kno’wn density) and X-ray flux was 
measured as part of the "Antarqui" campaign, ion density 
versus time profiles will be constructed and compared to 
observed ion densities. These are, in turn, proportional 
to conductivity enhancements. This will enable us to 
obtain a qualitative and quantitative grasp on, and there- 
fore a better understanding of, some of the processes and 
ionisation sources of the low latitude middle atmosphere. 



CHAPTER II 


DESIGN AND DISCUSSION OF INSTRUMENTS 

2 . 1 Ion Collection Theory 

Both types of instruments used in these studies 
employ a similar theory of ion collection; this simple 
theory o£ positive ion collection will now be briefly dis- 
cussed. 

The fundamental equation for positive ion col- 
lection to a probe surface, before the instrument becomes 
saturated is: 

J = N’*'eS + uN’^e VV-eDVN'^ C2.1) 

where N is the positive ion density, J is the current 
density; D is the coefficient of diffusion of positive 
ions; V is potential; e is electrostatic charge; and S is 
the convection velocity (Farrokh, 1975). 

The first term of the right hand side of equation 
C2.1) is a convection term, the second represents the ef- 
fect of conduction, and the third term is due to diffusion. 

If it is assumed that the air flow is parallel to the central 
axis of the instrument, the convection term can be neglected 
and equation (2.1) becomes: 

J = yN'^eW - eDVN"* (2.3) 

Also, with the added assumption that ?N has a 
value at most the same order of VV, the second term of 
equation (2.2) can be neglected (Farrokh, 19"5) and equation 
(2.2) becomes: 


J 


yN' eVV 


( 2 . 5 ) 
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Equation (2.3) is the fundamental relation for a 
simple mobility charged particle collection m'^chanism. 

Hoiilt (1965) and Sonin (1967) have demonstrated that the 
equation of the form (2.3) is applicable to a number of 
different type probe devices and geometries operating in a 
conductive, continuous medium. 

2 . 2 The Blunt Probe 

Th;i bipolar conductivity of a weakly ionized 
plasma, such as the lower ionosphere in the 30-80 KM 
altitude range, can be measured by a two electrode, simple 
mobility mechanism device known as a blunt probe. This 
instrument was originally designed and constructed by 
L. C. Hale, D. P. Hoult and associated students of The 
Ionosphere Research Laboratory in 1964. Since then, it 
has been used regularly and extensively by this group in 
studying the "normal" daytime and nighttime D-region. Also, 
the effects of events such as solar eclipses, "v/inter anomalies", 
and polar cap absorptions, upon the lower ionosphere have 
been investigated v;ith the blunt probe technique (Mitchell, 

1973) . 

Hoult (1965) first described the theory of col- 
lecting charged particle species by a subsonic blunt probe 
under certain specific conditions. These conditions are 
that; 1} the convective flow must be much greater than 
diffusion and "ion mobility" until very close to the col- 
lecting surface, 2) the space charge density is low enough 
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that it may be neglected in calculating the electric field 
in the probe, and 3) the mobility must be definable in a 
physical sense. 

If these conditions are satisfied, ion collec- 
tion can be expressed by a simple "mobility mechanism" as 
discussed above: 

dl = eNyEdA 

where dl is the current collected by the element of probe 
surface dA, E is the electric field at the surface, y is 
the ion mobility and N is the number density of charged 
particles collected. This equation is analogous to equation 
C2.5), the general equation describing simple mobility col- 
lection theory. 

Figure 2.1 is an example of the probe system. 

An aluminum disk is located at the end of the cylindrical 
payload housing. The probe disk faces downward, suspended 
from the parachute, into the flow. The central portion 
of the disk is used as the collecting electrode. It is 
held at the same potential as the rest of the disk, which 
now serves as a guard electrode. The remainder of the 
payload housing is used as a return electrode (Hale, 1965). 

The potential (V ) of the flat plate charged 
particle sensor, which consists of a central collecting 
surface and a concentric guard ring of radius a and R 
respectively, is varied. The electric field, E is given 
by : 


E 


C2Vp)/7tR 


C2-4) 
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Figure 2.1 The Blunt Probe Assembly 
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and the current reaching the collecting surface is: 


I + ^ 2a^V„a+/R 

p~ p - 


(2.S) 


where are the conductivities in the ambient plasma. 
Conductivities are determined from dlpj^/dV^, rather than 
Ip+/Vp, to remove drift errors in the instrument, is 
produced by applying a linear sweep voltage (V^3 between 
the return electrode and the probe. If is greater than 
a few volts. 




dV„ 


( 2 . 6 ) 


and hence, equation (2.5) becomes; 

R dl. 


0T 


2a' 


dV” 




(2.7) 


Calibration and conductivity data are converted 
from telemetry tape to a strip chart recording of current 
as a function of time. dIp_+/dV^ is computed from the slope 
of these curves and is given by: 


dl + 
„JP- 
dV, 


^ (dlpj/dt) 


( 2 . 8 ) 


3 - X3T7dt)^l 

where is the calibration resistance. This equation 

refers to a ground-based calibration in which the sweep 
voltage is applied through to the electrometer. 

Therefore, equation (2.7) becomes: 


C'+ 


R 


2a^R 



(2.9) 


cal 


cal 
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enabling us to obtain a+ independent of the sweep ampli- 
tude. 

Figure 2.2 shoves typical voltage and current 
characteristics of the probe experiment. 

Since radar constantly observes the probe's 
positionj the conductivity can be plotted as a function 
of altitude, Z. 

The main advantage of the blunt probe, over 
slender body geometries such as the Gerdien condenser is 
that its collection efficiency is much less sensitive to 
angle of attack on descent and, therefore, less likely to 
be affected if the probe package is swinging at the end 
of the parachute after deployment. 

The most significant problem with the blunt 
probe technique is that it measures only conductivity 
and not number density. A properly selected average 
mobility is a reasonable way to attack this problem (Mitchell, 
1973) . 

Another problem concerns contamination from the 
probe itself, i.e, by outgassed water vapor. The result 
would be to hydrated ions, thereby decreasing their 
apparent mobility and lowering the measured conductivity 
(Mitchell, 1973). However, ion conductivity measurements 
on geophysically similar days are similar, indicating that 
contamination due to water vapor is not significant (Mitchell, 
1973). 
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Finally, photo-emission on the collecting sur- 
face does not appear to be a significant source of error, 
since this surface is virtually always shielded from direct 
sunlight, 

2 . 3 The Gerdien Condenser 

The other simple type of instrument used in 
measuring the electrical properties of the middle atmo- 
sphere is the Gerdien condenser. This device, first 
constructed by H. Gerdien in 1905, consists of two con- 
centric cylinders positioned so that air flows through 
the apparatus. At sufficiently low potential difference 
between the inner and outer cylinder (Figure 2,3a), ions 
are collected at a rate depending upon the potential dif- 
ference across the cavity and flow rate through it. 

The theoretical current-voltage characteristics 
of the simple Gerdien condenser geometry as shown in 
Figure 2.3a are well knoiim. As in the blunt probe tech- 
nique discussed above, this device also uses a simple 
mobility mechanism. A sweep voltage is applied between 
the cylinders and an electrometer measures the resulting 
collector current. Figure 2.3b shoxtfs a typical I-V curve 
for a gas of single ion mobilities. 

If there are many type ions at different densi- 
ties and mobilities the I-V characteristics may be obtained 
by superposition. The current obtained before saturation 
is (Farrokh, 1975); 
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I = e§ V ! MnN„ 


( 2 . 10 ) 


where I is the current; e is the permeability of the 

2 TT El L( 

medium; n is the number of ion species; C = T h ' (^ 7l T ') 
the capacitance of the Gerdien condenser; L is the length 
of the cylinder; and Rq and are radii of the outer and 
inner cylinders. 

Since the number of ions in the flow sample 
is finite, the instrument has a saturation current, i.e. 
increased voltage will not produce a similar increase in 
current and the I-V characteristic becomes nonlinear. 

For a single type ion species, the negative saturation 
current I is proportional to IJ(N + N) where N is electron 
density and U is the average stream velocity. 

After saturation, equation (2,10) has the 

form: 


n 


= e 


( 2 . 11 ) 


- Vn I M N 
0 e 0 ^ n n 

where Vq is the saturation voltage. 

Figure 2.3c shows the I-V characteristics of 

a Gerdien condenser collecting tv/o types of positive ions. 

Here, the saturation current, Iq, is proportional to 

UCN^ + N 2 ) • The derivative of the current -voltage 

characteristic is proportional to ^ 2^2 ^0 ^ ^0. ' 

The ion conductivity, N~ e can be calculated 
from the slope of the I-V characteristic, since, from 
equation (2.10), ^ = a~ , and N— Lf can be obtained from 
the saturation current. N- or U must be known (or assumed) 
to obtain values of p+. 
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The Gerdien can be operated as a mobility 
spectrometer as demonstrated in Figure 2.3c. It is 
necessary to have an accurate I-V characteristic to 
distinguish ions o£ different mobilities. 

If ions of one polarity collect more quickly 
than others of opposite polarity, a space charge will 
be present between the two electrodes of the condenser, 
and may distort the electric field in the condenser. 
Farrokh (1975), however, shows that significant errors 
due to a space charge effect do not occur below an 
altitude of 80 Km, Photo-emission on the collector sur- 
face does not seem to present serious problems either. 

The obvious disadvantage of this device is its 
high sensitivity to the angle of attack, although better 
parachute design lessens the severity of this effect. 

The main point to consider, however, is that 
the Gerdien method gives a direct measurement of ion 
density and conductivity. The blunt probe technique 
yields the charged particle conductivities, but to deduce 
ion density, the mobility, u? must be knora. Since the 
ion mobilities in the lower ionosphere are unknown, the 
blunt probe technique is not reliable in obtaining densi- 
ty measurements. 
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2. 4 Comparison of Blunt Probe and Gerdien Condenser 

Measured cLn'ductivities 

As part o£ this research effort, blunt probe and 
Gerdien measured conductivities obtained during geophysi- 
cally similar days and times were compared. This was 
done to examine the reliability of not only the Gerdien 
positive conductivities, but also of the Gerdien positive 
ion density. The latter will play an important role in 
subsequent discussions presented in this report. 

Drawn on Figure 2.4 are the positive ion con- 
ductivity curves as measured by a blunt probe on 23 May 
at 1430 EST and a Gerdien on 28 May at 1630 EST. The ob- 
served values agree fairly well below 55 Km, but are 
noticeably different between 55 and 70 Bbn. Similar com- 
parisons of the Aladdin (1974) campaign data show much 
better correspondence between the blunt probe and Gerdien 
data (Leiden, 1976), The differences observed in the 
Peru data are still unexplained at this time. 

Fortunately, the measurement of ion number den- 
sity with the Gerdien is more reliable than conductivity 
measurements. This is true since, as discussed in Section 
2. 2-2, 3, the N values are derived from the saturation 
currents, while the positive conductivities are derived 
from the slope of the current-voltage characteristics 
which suffer from drift of zero errors. So, even though 
a difference in measured conductivities exist, the Gerdien 
number densities are probably more in line with the blunt 
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probe conductivity measurements and, therefore, it will 
be these positive ion densities that will be used in 
a theoretical analysis of the problem in a subsequent 


section . 



CHAPTER III 


EXPERIMENTAL RESULTS 

3. 1 Conductivity Profiles 

As part of the '*Antarqui” D-region study, held 
in May, 1975 in Peru, three blunt probe and one Gerdien 
condenser instruments designed and constructed here at 
The Ionosphere Research Laboratory under the guidance of 
Dr. L. C. Hale were flown in an effort to obtain conduc- 
tivity and density measuTe^.ents of the low latitude middle 
atmosphere. A large number of other measurements were per 
formed by other groups and regular meteorological observa- 
tions were made as part of this study, making it possible 
to obtain a "snapshot" of the upper atmosphere during chis 
period. This is especially important since these measured 
quantities are usually dependent upon a number of diverse 
parameters . 

The instruments were launched at different times 
at the Chilca Rocket Range, using Super-Areas sounding 
rockets as launch vehicles. The latitude and longitude of 
the experiment at the Chilca Rocket Range was 12°15’15"S 
and 76®47’20"W, respectively. The observed data is now 
presented and discussed in some detail. 

Figure 3.1 shows the measured conductivities 
for the three blunt probe measurements of the 28th (1430 
EST), 23rd (2007 EST] , and 24th (0200 ESI) of May respec- 
tively. (Appendix I lists the mini-computer program that 
calculates and plots these conductivities from the raw 
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data using the procedure outlined in Chapter II.) They 
correspond to a mid-day, late afternoon, evening and late 
night environments . Smooth curves have been constructed 
to best fit measured values. Only on the 1430 EST, 28 May 
flight was information obtained on negative conductivity. 
Figure 3.2 shows a comparison of the three blunt probe 
measurements with an average mid-latitude day (Hale, 1975] 
and a nighttime Ifhite Sands Missile Range positive con- 
ductivity profile (Mitchell, 1975). 


3. 2 Discussion of Data 

A number of significant points concerning these 
observations should be noted. On all the daytim:? profiles, 
a "knee" occurs at approximately 63 Km. This compares 
very well with man)' past measurements of positive con- 
ductivity (Hale, 1974) and probably represents physically 
the altitude at which Lyman « and cosmic ray ionization 
crossover; the cosmic rays being responsible for the pro- 
duction of ions in the lower, denser atmosphere and the 
Lyman a radiation ionizing NO above 63 Km. At night, 
production of ions due to direct Lyman a ceases and the 
conductivity curve "rises" upward indicating the loss of 
this ionization source. 

It will also be noted that the conductivity is 
approximately constant for the day measurements in the 
50-65 Km altitude range indicating that positive conductivity 
is relatively independent of atmospheric dep.sity. Cipriano 
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(1973), using statistical techniques, shows that the posi- 
tive conductivities correlate extremely well with atmospheric 
temperatures as measured by meteorological "rocket - sondes . ” 

He deduces the temperature coefficient of the conductivity 
increased from 1.7I/®K at 35 Km to 4.6%/°K at 58 Km. The 
positive conductivity variations do not correlate well 
with negative conductivity though (Hale, 1973), 

On the 0200 EST, 24 May profile, there appears 
to be wave-like structure in the data from 40 to 50 Km. 

These may represent equatorial gravity waves in the middle 
atmosphere: wind data and wavelength considerations bear 

this hypothesis out (Clark, 1976). Similar structures 
have not been observed at mid-latitudes, indicating that 
these short wavelength disturbances are an observable 
feature of the equatorial upper atmosphere, as predicted 
theoretically (Holton, 1975), 

The 1430 EST, 28 May positive and negative 
conductivity curves reveal an absence of electrons below 
40 Km where the profiles then diverge very rapidly, indi- 
cating the presence of many daytime electrons. 

The superposition of these profiles in Figure 
3.1 yields some interesting results also. The Peru con- 
ductivity curves are very similar up to 50 Km or so, then 
there is an enhancement of the late night curve. The day 
and evening curves remain fairly constant up to 65 Km 
where the difference in conductivities arises due to direct 
Lyman a being present in the day but absent at night. 


AWTARQUr 
PERU , !975 


1430 EST , 28 mV 

2007 EST » 23 MAY 



Figure 3.1 Conductivities for 1430 EST, 28 May; 

2007 HST, 23 ^fay; and 0200 EST, 24 May, Peru 
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(We shall see that although scattered Lyman ct is present 
during the night hours? it is insignificant as compared 
to other sources.) The night and evening curves inter- 
sect again, a bit above 75 Km, the region o£ maximum 
enhancement being at about 65 Km. 

This observation is important since the layer 
of enhancement, as will be shown in the next chapter, cor- 
responds very well with the altitude range of significant 
measured, ion production due to celestial X-rays. Here 
is evidence possibly substantiating the reported effect 
of celesticl X-ray sources upon VLF radio wave propagation 
(Sharma, et al, 1972). If this enhancement was, on the 
other hand, related to a return to nighttime steady state 
conditions, the layer affected is expected to be much 
deeper, perhaps down to 40 Km. This is certainly not 
^vhat is observed in Figure 3.1. 

Finally, two '^typical'' mid-latitude conductivity 
profiles (one day and one evening) were drawn for compari- 
son purposes, Figure 3.2. The "knee” still occurs at 
65 Km but the measured conductivities are greater by a 
factor of 2; this is due to increased cosmic ray production 
of ionization at more northerly latitudes, as predicted 
by theory. The observed nighttime enhancement is still 
about 1.5 times the day conductivities. Interestingly, 
the altitudes where the nighttime enhancement begins and 
ends is lower than those observed in the equatorial flights. 
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Positive Ion Densities 

Positive ion concentration was measured directly 
by a Gerdien condenser on 24 May, 1630 EST, as part of 
the overall campaign. This data is presented in Figure 3.3 
(Leiden, 1976), 

Since a comparison of the directly measured 
Gerdien ion densities and, the derived blunt probe densities 
is desired, an analysis of the derived positive ion con- 
centration from the blunt probe conductivity follows. If 
one assumes a small ion model, positive ion densities can 
be calculated from the measured conductivities by (Hale, 


1973) : 


N 


ep- 


(3.1) 


v/here e is the electronic charge and y the averaged posi' 
tive ion mobility- Now, the mobility of a single ionic 
species is (Dalgarno, 1962): 


p = 


TgP 


(3,2) 


where Pq and are the sea level reference pressure and 
temperature respectively, and p ^ is the reduced mobility, 
i.e. the mobility of the ion species at STP. P q assumed 
to be 1.8 X lO”^ m^ v"^ s"^ for positive ions (Loeb, 1955). 
Putting equation (3.1) into (3.2): 


N 


epoTPg 


C 5 . j ) 
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T and P are obtained from meteorological data 
or, in this case since the rocket-sonde data was ques- 
tionable, a standard atmosphere. The COSPAR International 
Reference Atmosphere as listed in the DNA Reaction Rate 
Handbook (1972) is used, making it possible to obtain 
positive ion density from only positive conductivity data. 

The N ’s were calculated as explained above for 
the 1430 EST blunt probe conductivity data; they are 
plotted on Figure 3.4. A number of features are to be 
noticed. The "knee" of the conductivity curve is re- 
flected in the ion density curve as a deep minimum. Again, 
this is interpreted to be the "crossover" region betv;een 
Lyman and cosmic ray production. In the 65-50 Km layer, 
the slope is approximately proportional to air density 
(CIRA, 1972) indicating large numbers of electrons with 
a density roughly constant in this range, satisfying the 
steady state lumped ion continuity equation: 

Q = + a.N'^N" 

del 

where Q is ion production rate and and are the dis- 
sociative and ion-ion recombination coefficients (Hale, 
1973). 

If one compares Figures 3.4 and 3.3, it is 
obvious that the Gerdien densities and blunt probe calcu- 
lated densities differ by about a half order of magnitude. 
It appears that the assumed reduced mobility, Uq, is in- 
correct since a small ion model does not seem to agree 
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with recent experimental results (Hale, 1976), and that 
obviously large ion models will now have to be considered. 

It is concluded then that the Gerdien N^’s are probably 
better estimates than those calculated by assumed mobili- 
ties, and, thusly, will be used in subsequent calculations 
requiring ion densities. 

3. 4 Experimental Errors 

In experiments of this type, discussing absolute 
errors in conductivity measurements is very difficult to 
say the least. The theory of positive ion collection is 
well understood and, therefore, a high level of confidence 
can be put upon positive conductivity data and subsequently 
derived parameters. Since the collection of negative 
ions and electrons on the blunt probe surface is a non- 
linear function of electric field/atmospheric pressure, 
negative conductivities are difficult to interpret in an 
absolute sense. 

Some general statements can be made about 
relative errors. Altitude accuracy is usually +1 Km or 
better for all data. Blunt probe data itself is accurate 
to +10"o for subsonic velocities. Current-slope measure- 
ments are accurate. Geometrical effects and errors 
in probe potential account for j,10% error. Therefore 
the absolute accuracy of the conductivity measurements is 
+ 25^0 (Hale, 1967), The reliability of conductivity ratios 
is much greater, probably within 10°5. Lastly, these errors 
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are very similar to those in the Gerdien experiment and 
the estimated accuracy is comparable. This indicates 
that the data reported does have real physical significance 
The most important fact to note, however, is 
that all of the curves draim or mentioned in this chapter 
have been obtained and reduced in consistent fashion, 
putting a much higher level of confidence on the values 
when used for comparative purposes, as we shall indeed do. 


CHAPTER IV 


THEORETICAL CONSIDERATIONS 

4 . 1 Ionization Sources 

The significant sources of lower D-region 
ionization were briefly mentioned in Section 1.1, In 
this section, the actual nighttime production rates xv’ill 
be discussed and their relative importance assessed. 

Figure 4.1 (adapted from Rowe, 1972) presents 
the ion-pair production rate of these various sources of 
nighttime middle atmosphere ionization. The curves repre- 
sent the production rates as a function of altitude. For 
the purposes of this research, all the sources will be 
assumed to be isotropic, although the galactic X-rays may 
be sensitive to the zenith angles of intense point sources 
and, possibly, the upper atmosphere itself may act as an 
X-ray telescope, making the X-radiation non- isotropic . 

At 65 Km, the cosmic ray production, which is 
moderated by the solar activity cycle but fairly constant 
over shorter time periods, and the balloon measured 
galactic X-ray production (Barcus, 1975} are equal and 
represent the most important knovm ionization sources at 
this level. Scattered Lyman a, another possible ionization 
source, would contribute less than 10^ of the production, 
as compared to the above mentioned sources, at 65 Km and 
drops off very steeply at lower heights. 
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The X-ray curve decreases rapidly and ceases to 
be a factor below 40 Km or so. The cosmic ray production 
increases, proportional to air density, to much lower 
altitudes (Goldberg, 1975; Velinov, 1968). Finally, 
the mid-latitude precipitating electron production rate 
indicates that this is certainly not a significant source 
below 70 Km (Potemra and Zmuda, 1970). 

From these considerat ions , it appears that the 
galactic X-ray source should have a definite influence 
upon middle atmosphere ionization, especially in the 70-50 
Km altitude range. However, the matter of an enhancement 
in ion production due to X-rays is still controversial 
since most studies h? ce investigated the effects on VLF 
phase by X-ray sources and this phenomena occurs at a much 
higher altitude, about SO Km wiiere ionization is due to 
a number of sources and it is therefore more difficult to 
distinguish among them. With this in mind, an attempt 
has been made to model, in a simple and st raightforvvard 
manner, the ionization processes in the lower D-regicn, 
where it is felt the X-ray source would have a greater 
effect on production than at 80 Km. The theoretical be- 
havior of the model can be studied, and then compared to 
the measured conductivity and positiv'e ion number density 
'Miurnal change” in the region of interest, using the ion 
densities of Figure 3.4 as an initial condition. 
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4. 2 The Continuity Equation 

In the lower ionosphere, transport of positive 
ions can be neglected since the lifetimes of the ions are 
small as compared to transport times of neutral wind 
velocities (Rishbeth and Garriot, 1969^. Therefore the 
lumped, effective positive ion continuity equation may 
be written: 


dN' 

dt 


a 

d e 


a^N‘N 


C4.1) 


where N is the positive ion density; Q is the total ion 
production rate; is the dissociative recombination 

coefficient; and is the ion-ion recombination coef- 

ficient. 


Under nighttime conditions there are no 
electrons present and equation (4.1) becomes: 


dN 

dt 


= Q - a.N'^N" 


Finally, assuming N‘ = N we have: 


dN 

3F 


Q “ a^N 


r 2 


(4.2) 


(4.3) 


This equation will now be solved and then used to cal- 
culate the time dependence of the positive ion concentra- 
tion. 


4. 3 


Solution of the Continuity Equation 
Equation (4.3) has the general solution 


N (t) 


1/2 

C^) tanh[tanh"^(nQ(^)) + 
“i ^ Q 

t(a.Q)^/^] 


(4.4) 
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T* 

where is the initial value of N (positive ion density) 
and here serves as a boundary condition of this differential 
equation (Nawrocki and Papa, 1961) . 

Equation (4.4) can be solved numerically by 
logarithmic expansions of the form: 

C4.5) 


tanh'^ X = -|-ln 


if -1 < X < 1, and 


tanh X ~ 


e - e 

e + e 


(4.6) 


4.4 


Analysis of Calculated ion Density Versus 
Time Profiles 


The N versus time profiles for each of the 
altitudes 55 - 65 Km were calculated using the equations 
of Section 4.3 and the Gerdien measured densities at 
1630 EST, 28 May as n^, the initial ion density. Tv;o 

values of Q, Q„smic ray ^cosmic rayrX-ray’ 
varying reasonable values of were employed in this 

procedure. Four parameters in these calculations are 

presumed to be well known, ^ » Qv the 

initial and final positive ion concentrations. Here, 

it is assumed that the measured X-ray production rates, 

and this measurement was made only once, pictured in 

Figure 4.1 are constant. Since the 1430 EST and 2007 EST 

conductivities had little variation below 65 Km, it was 

assumed that the 1650 EST Gerdien densities would be 

reliable initial conditions. The final ' s , however, 
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were obtained from the reasonable assumption that N "ucr 
(Hale, 1973) and therefore the increase of N was calcu- 
lable from the observed enhancement of positive con- 
ductivity between 1430 EST and 0200 EST. It is also 
assumed that by 0200 EST a steady state, or near steady 
state, condition prevailed in the height range under 
consideration. With these conditions and assumptions, 
is the only free parameter. 

Figures 4.2 - 4.12 show the result of these 
calculations. Each solution asymptotically approaches 
a steady state ion density, N given by: 


N 


ss 


C^) 

^ a . 


1/2 


(4.7) 


For each of the two values of Q, as labeled, 
at a given altitude there appear to be a range of that 
will make the solution approach the observed enhancement 
of positive ion density, also labeled, within the limits 
of experimental error. 

Further, there will be two pairs of Q and a 
that will force the solution to approach asymptotically 

'i' 

the final or near final N observed and have relaxation 

s s 

times also on the order observed, v8 to 10 hours. These 
’’best fit" values from each of the eleven altitudes con- 
sidered are plotted on Figure 4.13. 

Clearly, there is a large difference in the 
values of with and without the added to the 

1 A IT y 


'COS rav 


Below 58 Km, increases suddenly in both 


cases, this seeming to indicate a new chemistry is occur- 
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ring. The fairly large scatter in is due to the scatter 
in the Gerdien observed The suppression of from 

61 to 59 Km results from an enhancement of N"*" in this layer 
(Figure 3.3) which may be an effect due to the temperature 
dependence of ion density (Cipriano, 1973) or of some 
unhnow ionisation source [see below). We v/ill now pre- 
sent an argument for the latter explanation. 

It is important to note, however, that the 

of the cosmic ray + X-ray production calculations is in 

closer agreement to the ion-ion recombination coefficient 

as given by the DNA Handbook (1972). This value, 

- 7 3-1 

3 X 10 cm s IS almost equal to the average value of a. 

1 

-I ” ^ ^ 1 

in the 65-58 Km range, 5.1 x lO^'^cm'^s’ , The average 
value for the cosmic ray only production late calculations 
is below 2 X 10 'cm'^s"'^. 

It appears that the cosmic ray + X-ray produc- 
tion model is more consisf^nt with the accepted and 
standard values of two body ion chemistry reaction rates. 

This conclusion is not totally surprising 
because during this rocket campaign. X-rays were measured 
at these levels and found to be comparable to calculated 
cosmic rays as ionization sources. Further, the observa- 
tion made in Section 5.2, showing that the vertical ex- 
tent of nighttime enhancement is coincident v/ith the layer 
of maximum ion production due to X-rays, supports the 
conclus ion. 
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This stud/ indicates that the presence of 
X-rays should be considered as an important source of 
ionization in the lower D-region and that some of the 
observed radio wave effects attributed to transient celes- 
tial X-ray sources can be verified theoretically. 

4 . 5 A Constant Recombination Coefficient Assumption 

Previous measurements suggest that the ion- ion 
recombination coefficientj a^, is constant with height 
between 55-70 Km (Hale, 1976). Therefore, equation (4.4) 
was employed again; this time was held constant and Q 
was varied to obtain the proper, observed, steady state 
positive ion density. Figure 4.14 is a height profile of 

additional ion production, AQC = (a^N ) - Q^os ray + X-ray^ " 

that is required to explain the nighttime enhancements if 

one assumes a- to be a constant. Three values of a. were 
1 

assumed; 4, 5, and 6 x lO'^cm^s'^. It is evident from 

the figure that under this assumption, a large ('vl.5 x 

Q V sharply peaked ai: 60 Km, additional layer 

^cos ray + X-ray-^ j r / i/ 

of ionization production is needed to explain the "Antarqui” 
data. The possible physical significance of this calcula- 
tion remains unknown at this time and has not been ex- 
plored further in this research. 
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Figure 4.14 Additional Production Needed 
I£ Constant is Assumed 



CHAPTER V 


CONCLUSIONS 

The ionization of the low latitude middle 
atmosphere has been investigated by the coordinated D- 
region ’'Antarqui" experiment, 25-28 May, 1975, Peru, 

Using data obtained during this campaign by The Pennsylvania 
State University Ionosphere Research Laboratory rocket 
launched parachute borne blunt probe and Gerdien condenser 
conductivity measuring devices, the present study concludes 
that : 

Cl) The relative constancy of positive con- 
ductivity between 65-50 Km indicates that this parameter 
is not strongly dependent upon atmospheric density. This 
conclusion has been stated by others using mid-latitude 
conductivity data (Mitchell, 1973). 

(2) The "knee” or minimum of the day, evening, 
and night conductivity curves occurs at about 65 Km; this 
is the same altitude that the "knee" appears on mid- 
latitude data (Hale, 1973). 

(5) Low latitude conductivities (and ion densi- 
ties) are lower than those at mid-latitudes. 

(4) Positive ion densities calculated by assuming 
a reduced mobility and observed blunt probe conductivities 
indicate that in the 65-50 Km layer they are approximately 
proportional to air density. This indicates a large 
number of electrons with roughly constant density, satisfyin 


• r . 
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a steady state, lumped ion continuity equation (Hale, 1973). 

(5) The region o£ enhancement coincides very 
xvell with the altitudes where the ion production due to 
celestial X-ray sources are assumed to be on the order of 
other ionization sources. If this enhancement was merely 
a return to steady state nighttime equilibrium conditions, 
the affected layer would be signi f icantly thicker, perhaps 
down to 40 Km, than the 50-70 Km region of observed 
enhancement (Hale, 1976), 

(6) The observed day to night variations of 
positive conductivity indicated that, with conclusion (5) 
above, cosmic ray ionization below 65 Km is at least neces- 
sary to account for these variations. Further, theoretical 
calculations indicate that an additional source of ioniza- 
tion, due to X-rays, is significant (i.e. comparable to 
cosmic ray ionization) in the 65 ^5 Km height range, if 
one assumes an ion-ion recombination coefficient of 

3 X 10~'crn‘^s” . This conclusion possibly is related to 
some of the reported changes in phase of VLF radio waves 
by transient, celestial X-ray sources. These individual 
sources contribute significantly to ti total X-ray flux. 

(7) If one assumes (ion-ion recombination 
coefficient) is constant in the height range 55-70 Km, 

an additional, large, sharply peaked at 60 Km, ionization 
layer is required to explain the ''Antarqui" observations. 
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APPENDIX B 

PROGRAM TO CALCULATE N^ 
VERSUS TIME CURVES 
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